function [b,se] = persistence(x,n)
% n is the number of observation between the current variable and the
% future variable (i.e. 4 if quarterly data and we want annual
% persistence).

if nargin < 2
    n = 4;
end
[b,se] = beta_hac(x(n+1:end),x(1:end-n));
end

